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We reinvestigate the bandwidth-control and doping-control Mott transitions (BCMT and DCMT) 
from a spin liquid Mott insulator to a Fermi liquid metal based on the slave-rotor representation of 
the Hubbard model, Q where the Mott transitions are described by softening of bosonic collective 
excitations. We find that the nature of the insulating phase away from half filling is different from 
that of half filling in the respect that a charge density wave coexists with a topological order (spin 
liquid) away from half filling because the condensation of vortices generically breaks translational 
symmetry in the presence of "dual magnetic fields" resulting from hole doping while the topological 
order remains stable owing to gapless excitations near the Fermi surface. Performing a renormal- 
ization group analysis, we discuss the role of dissipative gauge fluctuations due to the Fermi surface 
in both the BCMT and the DCMT. 
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I. INTRODUCTION 

Landau- Ginzburg- Wilson (LGW) paradigm has been 
our unique theoretical framework for classical phase tran- 
sitions. Starting from an electron Hamiltonian, one can 
derive an effective LGW free energy functional in terms of 
order parameters associated with some symmetry break- 
ing. The LGW framework has been also applied to quan- 
tum phase transitions by taking temporal fluctuations of 
order parameters into account, usually called the Hertz- 
Millis theory. 

There are models of quantum phase transitions, on the 
other hand, which may defy interpretation in the LGW 
paradigm. Consider the superfluid-insulator transition 
of a boson Hubbard-type model. The boson density (n) 
and its phase (ip) are canonically conjugate, satisfying 
the uncertainty relation AnA(f> > 1, and the competing 
nature of the two variables results in the condensation of 
one variable or the other, depending on the ratio of phase 
stiffness and the compressibility. The quantum conjugate 
nature of the variables, satisfying an uncertainty rela- 
tion, lies at the heart of the quantum phase transition in 
this particular case. It is not obvious how LGW theory, 
written solely in terms of an order parameter, captures 
the inherent competing nature of the conjugate variables 
driving the quantum phase transition. 

As another example of a quantum phase transition 
where an order parameter description is likely to fail, 
we mention the metal to paramagnetic insulator tran- 
sition (Mott transition) found in the study of the two 
dimensional Hubbard model. As recent dynamical mean- 
field theory (DMFT) studies show, the transition is as- 
sociated with the vanishing of the spectral weight of the 
quasiparticle peak, but not with any symmetry breaking 
field, hence the order parameter approach in the LGW 
paradigm is not clear to be applicable. Q 

Recently, Florens and Georges (FG) reexamined a 
bandwidth-control Mott transition (BCMT) from a para- 
magnetic Mott insulator of a spin liquid to a correlated 
metal of a Fermi liquid at half filling in the Hubbard 



model, jj In order to describe the BCMT they intro- 
duced an elegant formulation based on the slave-rotor 
representation, and investigated properties of metallic 
and insulating phases of the model. Many of the prop- 
erties obtained at the mean-field level matched well with 
the more sophisticated DMFT calculations. (3|| In this for- 
mulation the competing nature of canonically quantum 
conjugate variables naturally appears. Within this the- 
oretical framework the Mott transition is understood by 
softening of bosonic collective excitations, physically as- 
sociated with zero-sound modes in a Fermi liquid. When 
these bosonic excitations are gapped, a paramagnetic 
Mott insulator with charge gap but no spin gap results, 
thus called a spin liquid. On the other hand, condensa- 
tion of the boson excitations causes a coherent quasipar- 
ticle peak at zero energy, resulting in a Fermi liquid in 
the low energy limit. 

In the present paper we investigate a doping-control 
Mott transition (DCMT) from the spin liquid to the 
Fermi liquid based on the slave-rotor representation of 
the Hubbard model. We find that the DCMT differs 
from the BCMT in the respect that the nature of the 
Mott insulator and the mechanism of the Mott transi- 
tion are different from each other. Hole doping results 
in a nontrivial Berry phase term to the boson field, lead- 
ing to an effective magnetic field for its vortex field in the 
dual formulation. It is shown that this effective magnetic 
field induces a crystalline phase of doped holes, coexisting 
with the spin liquid. On the other hand, the paramag- 
netic Mott insulator at half filling is the same spin liquid, 
but without any charge orders. We argue that the doped 
spin liquid with charge order evolves into the Fermi liquid 
via a continuous phase transition. 

The present scenario for the DCMT was discussed 
before, but based on the boson-only (Hubbard) model, 
where fermionic excitations are ignoredj^, 0| or decou- 
pled to the bosonic excitations in the renormalization 
group (RG) sense. In this paper we start from the elec- 
tron Hubbard model, and derive an effective bosonic field 
theory. It should be noted that this effective field theory 
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is totally different from that in the boson Hubbard model 
owing to the presence of damped gauge fluctuations re- 
sulting from gapless fermion excitations. Although the 
Berry phase plays the same role in both doped (boson and 
fermion) Mott insulators, nature of the Mott transitions 
would be different owing to the presence of dissipative 
gauge fluctuations in the slave-rotor representation of the 
electron Hubbard model. Performing an RG analysis, we 
show that the dissipative dynamics of gauge excitations 
makes both the BCMT and the DCMT in the electron 
Hubbard model differ from those in the boson Hubbard 
model. There also exists a previous study considering 
fermion excitations coupled to bosonic fields.^] However, 
this study starts from the quantum dimer model, and 
considers a valance bond solid instead of the spin liquid 
with a Fermi surface. Thus, the fermion excitations in 
the model are gapped (for the s — wave pairing case), 
thus ignored in the low energy limit. We would like to 
emphasize that our crystalline phase is nothing to do 
with Cooper pairs, 0,0,0] instead associated with doped 
holes. 



II. EFFECTIVE FIELD THEORY FOR THE 
MOTT TRANSITION 

We derive the slave-rotor representation of the Hub- 
bard model in the path-integral formulation. We note 
that FG derived it based on not only the canonical quan- 
tization method but also the path integral formulation. 
However, we argue that our path integral derivation 
more clearly shows the connection between Hubbard- 
Stratonovich (HS) flelds and rotor variables. 

We consider the Hubbard model in two dimensions 



H 



ija 



(1) 



Here f is a hopping integral of electrons, and U the 
strength of local interactions, iii = cl^Ciu is an elec- 
tron density. 

A usual methodology treating the Hubbard U term is 
a HS transformation. Using the coherent state represen- 
tation and performing the HS transformation, we obtain 
the partition function 

Z = J D[cicr,tpi]exp^- J dT(^c*„{dr - n)cicr 

:^cLc,,))j, (2) 



where is an order parameter associated with a charge 
density wave (CDW), and the chemical potential of 
electrons. Physically, the Lpi field corresponds to an ef- 
fective electric potential. In the usual mean-field man- 
ner the CDW order parameter is given by —i^i ~ 

Integrating over electronic excitations in Eq. (2) and 
expanding the resulting logarithmic term for the effective 



potential ipi, one can obtain an effective LGW free energy 
functional in terms of the CDW order parameter (pi . As 
mentioned in the introduction, it is not clear that this 
LGW theoretical framework has the competing nature 
of quantum conjugate variables because there exists only 
one CDW order parameter. One can say that the for- 
mulation Eq. (2) is exact, and thus the LGW framework 
may be a good starting point. However, an important 
point is how to expand the resulting logarithmic term. 
The expansion should be approximately performed, and 
thus one cannot say validityof the LGW framework for 
quantum phase transitions. |3| 

It is clear that the metal-insulator transition is as- 
sociated with charge fluctuations. One way controlling 
charge fluctuations is to introduce the canonical conju- 
gate variable of the charge density. Unfortunately, (pi 
is not the canonically conjugate variable of the charge 
density because it is an effective electric potential. 

We consider the gauge transformation for an electron 
field 



(3) 



Here e^'^' is assigned to be an annihilation operator of 
an electron charge, and fi^ an annihilation operator of an 
electron spin. In this paper we call e~**' and /^o- chargon 
and spinon, respectively. 

Inserting Eq. (3) into Eq. (2), we obtain 



Z 



dr 



J2 fUdr -fi- ldre,)h. -tYfl^^'~''^f^- 

-E(i^^?-^^^E/i/-))l- (4) 



Performing the HS transformation (1/4L/)(^| 



ULf + iLiLpi, and shifting tp^ into pi 



(y9j - drOr, we 



obtain the following expression for the partition function 

Z = j D[f.ia, e^, tpi,Li] exp|^- J dr 



ia ija 

+ Y}ULj -h iU{p, - drO,) - ip, Y 



(5) 



Integrating over the ipi field, one finds Li — J2a flafii^- 
ln this respect Li corresponds to the density variable of 
FG. 1] 

Eq. (5) has an interesting structure for the quantum 
phase transition. First of all, there is the competing na- 
ture of canonically conjugate quantum variables. The 
6i field is canonically conjugate to the charge density 
~ X^cr fiafiiy^ from the coupling term 

—iLidrOi of the Lagrangian derived above. These two op- 
erators satisfy the commutation relation [Oi^Lj] = iSij, 
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and thus the uncertainty relation AL^A^^ > 1 works. 
Fluctuations of the 9i field correspond to bosonic collec- 
tive excitations, here associated with zero sound modes 
of a Fermi liquid when it becomes condensed. 1] This can 
be justified from the fact that the dispersion of the di field 
in its condensed phase is given by that of sound waves. 

The quantity ipi is the CDW order parameter in Eq. 
(2). In the formulation presented in Eq. (5), however, 
it transforms as the time component of a U(l) gauge 
field. Under the U(l) gauge transformation for the mat- 
ter fields, ficr e^'^'fia and 9i 9i + the effective 
potential should be transformed into (pi ^ ipi -\- dr4>i- 
This gauge-field aspect of the order parameter is intro- 
duced due to the mapping of Eq. (3), which involved the 
new phase degree of freedom. 

Integrating over the Li field, Eq. (5) reads 

Z = I D[U,0,,^,]e-f''^\ 

i 

This expression is nothing but the slave-rotor representa- 
tion of the Hubbard model, obtained by EG in a different 
fashion. It is clear that the CDW order parameter ap- 
pears to be the time component of a U(l) gauge field. 
This can be understood by the fact that physics of the 
CDW order parameter is an effective potential. 

This effective Lagrangian should be considered to gen- 
eralize the LGW theoretical framework. If fiuctuations of 
the 6i fields are ignored, the resulting effective field the- 
ory belongs to the LGW framework. However, as clearly 
demonstrated by EG, 9i fluctuations are mainly respon- 
sible for the metal-insulator transition occurring in the 
Hubbard model at half-filling. Keeping the di fluctua- 
tions, the effective field theory for the Mott transition 
is naturally given by a gauge theory. In this respect 
the Mott transition should be viewed beyond the LGW 
paradigm. 

A standard treatment of the hopping term in Eq. (6) 
yields the effective Lagrangian 

icr 

~t ^ ^ {fiaPijfja + f jaPij fi<y) 
<ij><J 

i 

-t E ie^'-a^je-^'^ + e^^^ a^e"'^-), (7) 

where and f3ij are spinon and chargon hopping order 
parameters, respectively. 



A saddle point analysis results in the self-consistent 
equations 

{Y,flho) = l-5, (8) 

where 5 is hole concentration. 

Considering low energy fluctuations around this saddle 
point, one can set = ae*°'J , Pij = /?e"'J and Lpi = 
Tp, + where a = /i^/jo)! and /3 = |(e^«^e-«')l 

are amplitudes of the hopping order parameters, and ay 
and air are spatial and time components of U(l) gauge 
flelds. Inserting these into Eq. (7), we flnd an effective 
U(l) gauge theory for the Mott transition 

Lq = 2tNaf3, 

= E •^^o(^r - II -lip,- ia,r)fza 
<ij>cr 

i 

—2ta E cos(6'j — Oi — Qij), (9) 

<ij> 

where A'^ is a total number of lattice sites. Eq. (9) is our 
starting point for the metal- insulator transition. 

In this effective gauge theory two important facts 
should be taken into account since they discriminate the 
DCMT from the BCMT. One is an effective chemical 
potential /ie// — ^J' + Wi in the spinon Lagrangian L/. 
Particle-hole symmetry at half flUing causes the effective 
chemical potential to vanish. On the other hand, away 
from half filling the particle-hole symmetry is broken, re- 
sulting in a nonzero chemical potential. 

The other important feature is a Berry phase term 
arising from the phase-fiuctuation term in the chargon 
Lagrangian Lg 

= -E I dr[ — {dr9,)dr0,+l{Y,fU.)dreiy) 

At half filling the Berry phase does not play any roles 
because time reversal symmetry considered in this pa- 
per leads to (drOi) = 0, and the average occupation 

number of spinous is given by {J2a flaficr) — 1- Insert- 
ing this into the expression of Berry phase, one obtains 
Sb = -Y^i Jq dridrOt = -27ri^-gi, where is an in- 
teger representing an instanton number, here a vortex 
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charge. Thus, the contribution of Berry phase to the 
partition function is nothing because of e~'^^ = 1. Away 
from half filhng the Berry phase action is obtained to be 
Sb = — 27ri(5^j (ji with modular 27r. This results in a 
complex phase factor to the partition function, given by 



Zq^ where Q — qi is a total instanton 



number, and Zg, the partition function for a fixed Q. 
The observation of Berry phase gives the motivation for 
this paper. In this paper we investigate how the effect of 
Berry phase makes the DCMT differ from the BCMT. 



III. BANDWIDTH-CONTROL MOTT 
TRANSITION 

First, we discuss the BCMT. Zero effective chemical 
potential and no Berry phase effect result in the following 
effective field theory 



Le = ^ ^{drOi - airf - 2ta ^ cos(6'j - 9^ - )(! 



Recently, the present author examined deconfinement 
of fermions in the presence of a Fermi surface. It 
has been argued that the fermion Lagrangian Lr in Eq. 
(11) has a nontrivial charged fixed point [ij! as the 
quantum electrodynamics in (2 + l)D {QED3) without 
a Fermi surface. [1^, [l3 The present author investigated 
the stability of the charged critical point against instan- 
ton excitations.. 13j following the strategy in Ref. 16|. 
In the presence of a Fermi surface the conductivity a/ 
of fermions is shown to play the similar role as the fla- 
vor number N of Dirac fermions in the QED^.Jl^ Since 
the flavor number of Dirac fermions is proportional to 
screening channels for the gauge propagator, large fla- 
vors weaken gauge fluctuations in the QED3. In the 
same way the conductivity of fermions near the Fermi 
surface determines strength of gauge fluctuations. Re- 
markably, the charged fixed point is found to be stable 
against instanton excitations when the fermion conduc- 
tivity is sufficiently large. 13J This implies that the U(l) 
gauge field can be considered to be noncompact. Eq. (11) 
)ican be a stable theory against instanton excitations. In 
this respect the spin liquid state can survive beyond the 
mean field level. But, the spinous are not free particles 
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In the absence of U(l) gauge fluctuations this effective 
action was intensively studied by FG.jl] For the mean 
field treatment FG utilized large N generalization of the 
chargon field, and derived the saddle point equations in 
Eq. (8) at half filling. They found that there exists a 
critical U/t for chargon condensation. [Tlj In the case of 
U/t > (U /t)c chargons are gapped, but spinous are mass- 
less. Existence of charge gap but no spin gap corresponds 
to a spin liquid Mott insulator. In the spin liquid there is 
no coherent quasiparticle peak at zero energy, and only 
incoherent hump was found near the energy ±[/. In the 
case of U/t < {U/t)c condensation of chargons occurs, 
causing a coherent quasiparticle peak at zero energy in 
the presence of incoherent hump near the energy ±U. As 
a result a correlated paramagnetic metal appears. Fur- 
thermore, FG analyzed the saddle point equations near 
the Mott critical point {U/t)c, and obtained mean field 
critical exponents for the charge gap and the quasipar- 
ticle weight. They also found that the effective mass 
of quasiparticles does not diverge near the Mott critical 
point owing to the spinon dispersion. 

However, the mean field analysis of FG should be 
checked in the presence of U(l) gauge fluctuations since 
instanton excitations of U(l) gauge fields can cause con- 
finement of spinous and chargons, completely spoiling 
the mean field picture. In two space and one time di- 
mensions [(2 -I- 1)D] it is well known that static charged 
matter fields are always confined owing to instanton 
condensation. [T^ For the mean field picture of the spin 
liquid and the Mott transition to be physically meaning- 
ful beyond the mean field level, the stability of Eq. (11) 
should be guaranteed against instanton excitations in the 
RG sense. 



any more owing to long range gauge interactions, result- 
ing in an algebraic behavior of the spin-spin correlation 
function with an anomalous critical exponent. '15', '18^ 

The Mott transition beyond the mean field description 
is more complex owing to the dissipative nature of gauge 
fluctuations. Integrating over the spinous, one can obtain 
the effective action for the chargon and gauge fields in the 
continuum limit 



S. 



drd^r 



[drO - arf - 2tocos(V0 - a) 



+ J d-qr-^a^{qr,iuJn)Dj;l{qr,iuJn)a^{-qr,~iuJn), 

(12) 

where Di^i^{qr,iujn) is the renormalized gauge propagator, 
given by 



,y{qr,iuJn) = (Sf,i, - ^^^D{qr,iLL!n), 



D ^{qr,iujn) = Di^'^{qr,iujn) + 'n{qr,iuJn)- (13) 

Here Z)q ^(g^, iw„) = (g^ +u)^)/g'^ is the bare propaga- 
tor of the gauge field given by the Maxwell gauge action, 
resulting from integration of high energy fiuctuations of 
spinous and chargons. g is an internal gauge charge of 
the spinon and chargon. Il{qr,iujn) is the self-energy of 
the gauge field, given by a correlation function of spinon 
charge (number) currents. Since the current-current 
correlation function is calculated in the noninteracting 
fermion ensemble, its structure is well known fl^ l20l | 



Il{qr,iujn) = a{qr)\uJn\ 



(14) 



Here the spinon conductivity a{qr) is given by (j{qr) ~ 
ko/qr in the clean limit while it is (T{qr) ~ (Tq = kgl in 
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the dirty limit, where fco is of order kp (Fermi momen- 
tum), and I the spinon mean free path determined by 
disorder scattering. The diamagnetic susceptibihty x is 
given by X ~ where ruf ^ {tf3)~^ is the band mass 

of spinons. The frequency part of the kernel Il{q,iu!n) 
shows the dissipative propagation of the gauge field ow- 
ing to particle-hole excitations of spinons near the Fermi 
surface. 

Eq. (12) should be a starting point for the BCMT. 
In the study of FGQ U(1) gauge fluctuations are ig- 
nored, and thus the physical picture of the Mott tran- 
sition should be modified. In the absence of U(l) gauge 
fluctuations the transition falls into the XY universality 
class. However, long range gauge interactions alter the 
XY universality nature into the inverted- XY (IXY) uni- 
versality class if the Landau damping term in Eq. (14) 
is ignored, and only the Maxwell kinetic energy of the 
gauge field is taken into account. This means that if 
one considers a critical exponent v associated with the 
charge gap Ag ~ |J7 — Uc\'^ with the critical value Uc, the 
critical exponent changes from vxy of the XY transition 
to lyjxY of the IXY transition. Damped gauge interac- 
tions are expected to modify the IXY Mott transition. [2^ 

Performing the duality transformation for the phase 
field in Eq. (12), we obtain the dual vortex action 



a^, we obtain the effective vortex action 



Sv = drd r 



\id,-ic,)<t>\'+mlm' + ^\<f\' 



+U{d X c)l + —{d X c)l - iaf,{d x c)^ 
+ J dqr-af,(qr,iuJn)D~Xqr,iuJn)a^{-qr, 



(15) 



Here $ is a vortex field, and a vortex gauge field. m„ 
is a vortex mass, given by ~ {U /t)c — U /t, and Uy a 
phenomenologically introduced parameter for local inter- 
actions between vortices. {U/t)c is the critical strength 
of local interactions, associated with the Mott transition 
in the mean field level. 

In the dual vortex formulation the BCMT arises from 
controlling the vortex mass as a function of the parameter 
U/t. In the case of ml < {U/t > {U/t)c) condensation 
of vortices occurs, resulting in a Mott insulator of char- 
gons. In the case of ml > {U/t < [U/t)c) vortices 
are gapped, implying condensation of chargons, and a 
paramagnetic metal results. 

Performing the Gaussian integration for the gauge field 



Sy = / drd'^r 



\i^,-^c,)<^>\^ + ml\'^>\' + ^\'^>\ 



+Uidxc)l + —{dxc)t 



+ / drdTid rd ri-c^{r,T)Kf_„y{r - n^T - Ti)c^{ri,Ti), 



(16) 



where the renormalized gauge propagator Kf^^{r — ri, r- 
Ti) is given by in energy- momentum space 



K{qr,iujn) 



{q^^+uD/g' + a{qr)\uj„ 



(17) 



where g is a redefined variable including the susceptibil- 
ity. In the following we consider dirty cases characterized 
by cr(g,.) = (To- 

Before we analyze Eq. (16) by using an RG method, we 
consider two physical limits; one is erg —> correspond- 
ing to an insulator of spinons, and the other, ctq — > oo 
identified with a perfect metal of spinons. In the spinon 
insulator the kernel K(qr, itOn) becomes a constant value, 
making vortex gauge fluctuations (c^) gapped, thus ig- 
nored in the low energy limit. This is because long range 
gauge interactions (a^) make it massive the low energy 
mode (Goldstone mode) represented by the vortex gauge 
field, appearing at high energies. The usual action 
for the vortex field is obtained. On the other hand, in 
the perfect spinon metal gauge fluctuations are com- 
pletely screened by spinon excitations, causing the kernel 
to vanish, and the Maxwell gauge action for the vortex 
gauge field results. The resulting vortex action is reduced 
to the standard scalar QED3. Varying the spinon con- 
ductivity (Joj these two limits would be connected. 

We perform an RG analysis for Eq. (16). Anisotropy 
in the Maxwell gauge action for the vortex gauge field is 
assumed to be irrelevant, and only the isotropic Maxwell 
gauge action is considered by replacing [/, 1 /Ata with 
l/(2e^). Here is a vortex charge. In the limit of small 
anisotropy the anisotropy was shown to be irrelevant at 
one loop level. 23j To address the quantum critical be- 
havior at the Mott transition, we introduce the scaling 



e'r' and r 



and consider the renormalized 



theory at the transition point = 



Sy = y 

+ ^{d' X 



Z^\{d'^-zeyC,)<f\' + Z^^m^ 

(18) 
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where and are the renormaHzation factors de- 

fined by 

$ = e 2 c^ = e 2 'Zc'c^,., 

el = e-(^-^''Z-ie^,, = er^'-''^' Z^Z^^^m 

In the renormaUzed action Eq. (18) the subscript r im- 
plying " renormalized" is omitted for simple notation. 

Evaluating the renormalization factors at one loop 
level, the RG equations are obtained to be 



{A-D)el-(\N, + ^)et 
\ an/ 



dl 

duy 



(To - 

= (4 - D)uy + h{aa, el)eluy 



(20) 



Here A, C, P are positive numerical constants, and 
/i(cto, e^), g(CTo, e^) are analytic and monotonically in- 
creasing functions of ctq- is the flavor number of the 
vortex field, here given by Ny — 1. 

The first RG equation for the vortex charge can be 
understood in the following way. Integrating out critical 
vortex fiuctuations near the critical point ~ 0, we 
obtain the singular contribution for the effective gauge 
action 

5c = - ^ j d'^qr-Cf,{qr,iuJn)^f,i^{qr,i(^n)Ci,i-q,-iuJn) 



N„ 



qf + ^l 



-^^/q^+UJl + 



ao\UJn 



The first term in the kernel S((7r,ia;„) results from the 
screening effect of the vortex charge via vortex polariza- 
tion, causing the —XNye^ term in the RG equation while 
the second originates from that via spinon excitations, 
yielding the — (C/cro)e^ term. The first (4— D)e1 term 
denotes the bare scaling dimension of the vortex charge. 

For the second RG equation, unfortunately, we do not 
know the exact functional forms of h{ao,el) and g{cro, el) 
owing to the complexity of the gauge kernel. Owing to 
the spinon contribution K{qr,iu!n) [Eq. (17)] the kernel 
of the gauge propagator (c^) 



Dc{qr,iuJn 



1 



ql + ujI + elK{qr,iUn) 



(9r + Wn)(e? +cro|w„|) 



should be utilized instead of the Maxwell prop agator in 
calculating one loop diagrams. pH l23l l25l l26l| Note the 
dependence of the vortex charge in the effective gauge 
propagator. This gives the dependence of the vortex 



charge to the analytic functions /i(cro,e^) and g(cro,e^). 
Although the exact functional forms are not known, the 
limiting values of these functions are clearly revealed. In 
the limit of cto — the gauge kernel vanishes, thus caus- 
ing h{(jQ 0, e^) — > and g{ao — + 0, e^) 0. In the 
small (To limit the gauge kernel is given by 



Dc{qr,iuJn) 



Co 



6^ q^ 



thus resulting in /i((To,e^) — ChOa/el and .g((To,e^) — 
Cgo'^/e'^, where and Cg arc positive numerical con- 
stants. On the other hand, in the limit of cto ^ oo the 
gauge kernel is reduced to the Maxwell one Dc{qr, i^n) ~ 
l/(9r + ^n)- Thus, /i(cto ^ oo, 6^) Ci and g((To 
oo, e^) 02 are obtained, where ci and C2 are positive 
numerical constants. [2ll \2a^ As a result, Eq. (20 ) is re- 
duced to the RG equation of the <I>^ theory |27| in the 
limit of (To ^ 



duy 

~dr 



(4 - D)uy - p{Ny + A)u 



and that of the scalar QED3 [21|, [23 in the limit of (Tq 

00, 

''-^ = [A-D)el-\Nyel 

^ = (4 - D)uy + cieluy - p{Ny + 4)u^ - cze^ 
In the small (To limit the RG equations (20) result in 



del 
dl 
duv 

~dr 



{^^D)el~{\Ny + ^) 

V (Tq/ 



= (4 - D)uv -f Ch<JoUy - p{Ny + A)ul - Cgal 



In the scalar QED3 there is a delicate issue about the 
existence of the charged fixed point (e*^ ^ 0).[^ l^^, 
i26l | In this paper we do not touch this issue. Instead 
we assume the existence of the charged critical point in 
the scalar QED3 by controlling the A value. Then, the 
charged critical point (e*^(cro), u*[e*2(cro)]) in Eq. (20) is 
expected to vary as a function of the spinon conductivity 
in the range of 



e,f ((To 0) 

* r *2 / 



< el^ao) < el^aa 

o)l < u: 



f ((To)] <<[ef ((To 



(x.) = l/(AiV„), 
0)], 



where the fixed point (e*^(CTo ~^ 0), u*[e*^(CTo 0)]) cor- 
responds to the IXY one in the original boson model [Eq. 
(12)], and the fixed point (e*^(cro 00) , u*,[el'^ {ao 
00)]) coincides with the XY one in Eq. (12). The spinon 
contribution (cto) connects the XY fixed point to the IXY 
one smoothly in the chargon action Eq. (12). (2^ This im- 
plies that the critical exponents near the Mott transition 
change continuously, depending on the spinon conductiv- 
ity. This would be measured in some experiments. Be- 
cause the spinon conductivity depends on disorder, we 
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would have some different critical points by controlling 
density of disorder, resulting in various critical exponents 
between the exponents of the XY and IXY transitions. 
However, one interesting possibility should be taken into 
account that the glassy behavior of the chargon field can 
originate from random potentials. This important sub- 
ject is under current investigation. 

IV. DOPING-CONTROL MOTT TRANSITION 

Next, we investigate the DCMT, described by the ef- 
fective field theory 

ia 

-t(3 J2 (/*.e-"-/,, + c.c.), 

<ij>a 

= ^ '^{drOi - airY - "^ta ^ cos(6'j - 9i - aij) 

i <ij> 

+i5{dre, - a„). (21) 

Note the presence of the effective chemical potential and 
Berry phase. This Lagrangian is analyzed by employing a 
duality transformation. ,6j In the dual formulation the 
effect of Berry phase is represented as effective magnetic 
fields for dual vortex variables. 

Following the previous section, the duality transforma- 
tion of the chargon Lagrangian results in 

C,, = \{d^- ic^m^ + ml\^^ + ^|<i>|4 - h{d X c\ 

+[/(axc)2 + J-(axc)^, (22) 

where the U(l) gauge field was ignored in the mean 
field level. The Berry phase effect is reflected as an effec- 
tive magnetic field h = —2US for the vortex field in the 
term —h{d x c)r- Remember the expression of the vortex 
mass m.y ^ {U /t)c — U/t. A cautious reader may suspect 
that the vortex mass should depend on hole concentra- 
tion. From the discussion below Eq. (9) it is important 
to note that the effect of hole doping appears only in the 
chemical potential and Berry phase terms. Furthermore, 
at half filling Eq. (22) should be dual to the chargon La- 
grangian Lg in Eq. (11). Thus, the vortex mass should 
depend on only the parameter U/t. 

In the dual vortex formulation the BCMT is driven 
by controlling the vortex mass, as shown in the previous 
section. On the other hand, the DCMT is nothing to do 
with the vortex mass. Instead, controlling the effective 
magnetic field causes the Mott transition. This leads us 
to consider that the nature of the DCMT differs from the 
BCMT. 

The presence of the effective magnetic field reminds us 
of a well known Hofstadter problem for vortex fields. In 
the context of a superfluid-insulator transition this was 
extensively studied in Re fs. |5|,|Q]. Here we briefly sketch 



the procedure and key results. We first investigate the 
nature of a doped Mott insulating state in a mean field 
fashion, i.e., the absence of U(l) gauge fluctuations a^, 
and discuss the DCMT beyond the mean fleld level. 

Following Ref. 'B'l, we consider commensurate hole 
concentration 5 = p/q, where p and q are relatively prime 
integers. Under this effective magnetic fleld (5, the vor- 
tex Lagrangian Eq. (22) has g-fold degenerate minima 
in the magnetic Brillouin zone. Low energy fluctuations 
near the g-fold degenerate vacua are assigned to be '^i 
with Z = 0, — 1. A key question is how to construct 
a LGW free energy functional in terms of the 5*; fields. 
Constraints for an effective potential of 4'; are symmetry 
properties associated with lattice translations and rota- 
tions in the presence of the effective magnetic fields. 
Based on the symmetry properties one can construct a 
LGW free energy functional of 4";, and perform a stan- 
dard mean field analysis. In this free energy a superfluid 
of original bosons (chargons) is given by (^E";) — for all 
I — 0,...,(7 — 1 while a Mott insulator is characterized 
by (^;) 7^ for at least one /. Although the free en- 
ergy functional has all symmetries, the ground state can 
be symmetry-broken. In other words, the Mott insulator 
can have broken translational symmetries. 

To see this, one can construct a density wave order pa- 
rameter by considering bilinear and gauge-invariant com- 
binations of the low energy vortices . Condensation of 
leads to a nonzero value of the density wave order pa- 
rameter, causing a vortex density wave. A density wave 
of vortices can be interpreted as a crystalline phase of 
doped holes in the original language. In appendix 

we review the simple g = 2 case. Combining this chargon 
physics with the spinon physics, we can conclude that a 
doped Mott insulator consists of a density wave of char- 
gons and a spin liquid of spinous with a Fermi surface. It 
should be noted that this doped Mott insulator is differ- 
ent from the Mott insulator at half filling because there 
is no charge order in the undoped Mott insulator. 

Remember that the crystalline phase of doped holes 
is nothing to do with the spin liquid in the mean field 
level. Integrating out the gapped chargon degrees of 
freedom, we obtain the same spinon-gauge action for 
the doped spin liquid with that for the undoped one 
beyond the mean field level. It has been argued that 
the spinon-gauge action is a critical field theory at the 
nontrivial charged fixed point .[Tsll as discussed in the 
previous section. Thus, the spin-spin correlation func- 
tion shows a power law behavior with respect to fre- 
quency and temperature, Il8|| On the other hand, no 
infrared response for charge fiuctuations is expected ow- 
ing to the Mott gap. Instead, the charge order would be 
reflected in the electron density of states as a spatially 
modulated pattern because of the translational symme- 
try breaking. |2^ The electron density of states is pro- 
portional to the imaginary part of the electron green 
function, given by convolution of the spinon and char- 
gon propagators in the slave-rotor formulation. Thus, 
the spatial inhomogeneity of the chargon distribution re- 
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suits in the spatially modulated pattern in the density 
of states. Because there is excitation gap in the chargon 
spectrum, only incoherent hump would be shown in the 
electron spectral function. This is another different point 
from the usual density wave. 

One cautious person may suspect the coexistence of 
the spin liquid and charge density wave (CDW) because 
such a commensurate CDW can destroy the spinon Fermi 
surface through a space-dependent effective chemical po- 
tential, causing the spin liquid to be unstable. However, 
we argue that the spinon Fermi surface can be preserved 
even in the commensurate CDW when the Fermi surface 
nesting is not perfect due to interaction or frustration ef- 
fects. Considering the low energy vortex excitations 
near the g-fold degenerate vacua, one can find the effec- 
tive dual action 



H-iiPi- iair)fi, 



-tp J2 (fl^^'-'foo + h.c.)] cos(5 X a) 



<ij>cr 



<nm>/ 



(23) 



<tii'> 



Here I = 1, q corresponds to a color index of low energy 
vortex fields, and Fd^n"*!) is an effective vortex poten- 
tial determined by symmetry properties, where the coef- 
ficients are effectively doping dependent (see appendix). 
The last gauge action in the spinon sector originates from 
high energy contributions of matter fields, where g is an 
internal gauge charge of spinous. 

The question is what happens in the Fermi surface 
when vortex condensation occurs, resulting in transla- 
tional symmetry breaking. Ignoring spinon-gauge fluctu- 
ations a^i, as the mean field approximation, the spinon- 
gauge action is completely decoupled from the vortex- 
gauge action, as discussed before. This indicates that 
the Fermi surface is not affected by the CDW formation 
in the vortex sector. Now, we allow spinon-gauge exci- 
tations. Integrating out air in the limit of — > c», one 
obtains the constraint 



(V X C)^ 



When the vortices are condensed to cause translational 
symmetry breaking, the above quantity should depend 
on positions. This effect can be introduced in the spinon 
action by allowing a position-dependent effective chem- 
ical potential, interpreted as a higher order effect due 
to gauge fluctuations. In this case the commensurate 
CDW can destroy the Fermi surface. However, it should 
be noted that this depends on the shape of the Fermi 
surface. When the perfect nesting of the Fermi surface 
does not appear due to interaction or frustration effects. 



only partial parts of the Fermi surface would open the 
CDW gap, and other parts of the Fermi surface, not con- 
nected by the CDW wave vector, are expected to remain 
gapless. This would indeed happen when there is frus- 
tration, destroying the Fermi surface nesting. This ex- 
pectation coincides with our ignorance of spin ordering 
because strong frustration kills magnetic ordering. 

We should emphasize that the above discussion is ap- 
plied to the g ^ oo limit. Since we are considering low 
energy fluctuations, high energy matter flelds should be 
integrated out, resulting in the Maxwell gauge action, 
the last term in Sf of Eq. (23). Then, the above con- 
straint cannot be used because one cannot integrate out 
ttir directly. In this case of spin-charge separation due to 
a finite value of g, the position-dependent dual magnetic 
flux is not directly related with the spinon density owing 
to the spinon-gauge flux. If one utilizes average values in 
the constraint equation instead of the operator identity, 
the same argument above can be applied. Generically, 
there should be a gapless Fermi surface, not connected 
by the CDW wave vector, at least in the frustrated lat- 
tice. 

The present doped spin liquid is an interesting new 
phase in the respect that a conventional order described 
by the CDW order parameter (in the LGW paradigm) 
and an exotic order associated with a conserved internal 
gauge fluxfsol l3l| coexist. This phase is expected to be 
stable beyond the mean fleld level because instanton ex- 
citations can be suppressed via spinon excitations near 
the Fermi surface, as discussed before. A next important 
question is which phase this doped spin liquid evolves 
into. More concretely, when chargon condensation oc- 
curs, does the CDW order survive? Remember that in 
the BCMT chargon condensation results in a Fermi liq- 
uid, where condensed chargons are conflned with spinous 
to form electrons (quasiparticles). This corresponds to a 
Higgs-confinement phase in the context of gauge theory, 
where the internal gauge flux is not conserved. |30l Isij 
In the DCMT the chargon condensation also causes elec- 
tronic quasiparticles. In this case the CDW order is ex- 
pected to disappear, thus resulting in the same Fermi 
liquid as that in the BCMT. See Eq. (22). Because 
we are considering gapped vortex excitations, they can 
be ignored in the low energy limit. Thus, there remain 
uniform effective magnetic fields. This implies that con- 
densed chargons are homogenously distributed. As a re- 
sult the CDW order of chargons disappears. Another way 
to say this is that since there are no vortex charges (owing 
to the gap in the vortex excitations) in the Berry phase 
termEq. (10), the effect of Berry phase disappears. This 
is analogous to the case in the nonlinear a model for the 
quantum antiferromagnet, where the Berry phase effect 
can be ignored in the antiferromagnetic phase. [so) 

Now we discuss a critical field theory for the DCMT. 
Integrating over the spinon and gauge (a^) excitations in 
Eq. (23) as performed at half filling, one can construct 



the following effective field theory 

'' 1=0 

^(h - h,){d X c). + U{d X c)l + ^(5 X c)l 

+ J dTdTi(fr(fri^Cf_,{r,T)K^„{r - ri,T - ri)cy(ri, n). 

(24) 

Kf^i, results from the anomalous contribution of spinon- 
gauge fluctuations to vortex-gauge excitations, given by 
Eq. (17) with a different ctq owing to the chemical po- 
tential fj.eff- h = —2US is an applied effective magnetic 
field, and hq — —2U5q a nearby one with commensurate 
hole concentration 6q = p/q. One can estimate a criti- 
cal effective magnetic field he with a given U/t > {U/t)c 
by calculating the condensation energy. The critical hole 
concentration 5c corresponding to the critical magnetic 
field he would be different from 5q generally. In this case 
one may determine a moderate value of q near the critical 
doping 5c- Then, there remain residual effective magnetic 
fields h — hq, corresponding to the incommensurability 

5-5q. 

We propose that Eq. (24) is a starting point for the 
DCMT. If the vortex "superconductor" falls into the 
type-I class, the residual magnetic field would be expelled 
owing to the dual "Meissner" effect. A critical field the- 
ory for this Mott transition is expected to be without the 
residual magnetic field 

Scff^ / dTd\[Y^\id^-tc^)^SJi\^ + Vi^i) 

•' 1=0 

+C/(9xc)^ + J-(axc)2" 

+ j dTdTid^rd^ri]^c^{r,T)K^^{r - ri,T - ri)cy(ri, ti). 

(25) 

Because the effective vortex action depends on q and 
V{^i), it is difficult to predict critical vortex dynamics 
for general q values. The q = 1 case corresponds to the 
undoped spin liquid, already discussed in the previous 
section. In the q = 2 case the effective vortex potential 
is obtained to be 

Vii^i) = m'dV'ol' + iV-il') + "4(|V^o|' + l^il')' 
+Vi\^i\''\iP2f - vsUrM^ + H.c], 

well discussed in appendix. At the critical point = 
the last eighth-order term is certainly irrelevant owing 
to its high order. Furthermore, the cubic anisotropy 
term (^4) is well known to be irrelevant in the case of 
q < qc = 4. [23 As a result, the Heisenberg fixed point 
{vl = and m| ^ 0) appears in the absence of vor- 
tex gauge fluctuations. I23I Introducing the vortex gauge 
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FIG. 1: A schematic phase diagram in the slave-rotor repre- 
sentation of the Hubbard model 

fields at the Heisenberg fixed point, we have qualitatively 
the same fixed point with Eq. (16) except the q — 2 
vortices. Since the charged critical point depends on the 
spinon conductivity, the critical exponents vary as a func- 
tion of the spinon conductivity. At higher q values we do 
not understand the nature of the Mott transition ow- 
ing to the complexity of the vortex potential. Generally 
speaking, a continuous Mott transition from the U(l) 
spin liquid with a commensurate density wave order to 
the Fermi liquid is possible. 

On the other hand, if the vortex superconductor be- 
longs to the type-II class, the residual magnetic field can 
penetrate the superconductor, forming a dual Abrikosov 
"vortex" lattice. This corresponds to an incommensurate 
Mott insulator, where hole density is J ^ (5g.j^ In this 
case the nature of the Mott transition from the U(l) spin 
liquid with an incommensurate density wave order to the 
Fermi liquid is not clear owing to the Berry phase effect. 
Furthermore, the Landau damping term should be taken 
into account in the critical field theory as the case of the 
BCMT because it changes the nature of the Mott transi- 
tion. A continuous transition to the Fermi liquid may be 
possible in this case. A detailed analysis of this DCMT 
is beyond the scope of this paper. 

We propose a phase diagram Fig. 1 in the slave-rotor 
representation of the Hubbard model on two dimensional 
square lattice. Here UlSL, FL, and CDW represent U(l) 
spin liquid, Fermi liquid, and charge density wave, re- 
spectively. 

The route 1 is the BCMT from UlSL to FL at half 
fiUing while the route 2, that from UlSL + CDW to FL 
at commensurate hole concentration. In these cases we 
showed critical field theories, and discussed nature of the 
continuous phase transitions. 

The route 3 is the DCMT at a finite U/t. In this 
case the critical field theory depends on the nature of 
the vortex superconductor. Nature of the quantum phase 
transition from the UlSL with an incommensurate CDW 
to the FL is not clear owing to incommensurability. 
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V. DISCUSSION AND SUMMARY 



So far, we considered a zero flux state, and thus ob- 
tained the U(l) spin hquid with a Fermi surface. The uni- 
form spin liquid phase turns out to be stable (with respect 
to the flux phase) in the triangular lattice at the mean- 
field level near the undoped Mott transition point 0. 
However, it is important to consider a tt flux phase since 
this phase is usually obtained as a stable mean fleld state 
in the square lattice without frustration[32j . In the tt flux 
phase low energy spinon excitations are given by Dirac 
fermions near four nodal points. As a result the effec- 
tive spinon Lagrangian is obtained to be QED^. For 
the BCMT at half filling there is no dissipation in gauge 
fluctuations. The role of massless Dirac fermions is to 
weaken gauge interactions, resulting from screening of 
gauge charges owing to particle-hole polarization. Thus, 
increasing the flavor number N of Dirac fermions in the 
1/N approximation, the IXY transition is expected to 
turn into the XY one|3^. On the other hand, for the 
DCMT there remains dissipation in gauge fluctuations 
owing to a nonzero effective chemical potential. Thus, 
damped gauge fluctuations would still play some special 
roles in the DCMT. 

In this paper we discussed how the doping-control Mott 
transition differs from the bandwidth-control one based 
on the slave-rotor representation of the Hubbard model. 
We found that the doped Mott insulator consists of a 
crystalline phase of doped holes and a U(l) spin liquid 
with a Fermi surface while the Mott insulator at half flU- 
ing is the U(l) spin liquid without any charge orders. 
This originates from the fact that hole doping causes a 
Berry phase term to the chargon field. This Berry phase 
effect results in an effective magnetic field to the vortex 
field of the chargon field. In the dual vortex formulation 
we showed that the bandwidth-control Mott transition 
is driven by the sign change of the vortex mass while 
the doping-control one is achieved by the control of the 
effective magnetic field. The presence of dual effective 
magnetic fields leads to translational symmetry break- 
ing when vortices are condensed. We argued that this 
charge order does not destroy the spinon Fermi surface 
when there is strong frustration, causing the Fermi nest- 
ing to disappear. As a result, the spin liquid phase can 
remain stable to coexist with the density wave. Further- 
more, we pointed out that damped U(l) gauge fiuctua- 
tions resulting from spinon excitations should be taken 
into account for both Mott transitions because the na- 
ture of the Mott transitions is modified by the dissipative 
gauge excitations. Performing a renormalization group 
analysis, we showed that the Mott critical point depends 
on the spinon conductivity characterizing the strength of 
dissipation. This interesting result leads us to predict 
that varying the density of disorder would cause differ- 
ent critical exponents because disorder determines the 
conductivity of spinous. 
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APPENDIX: MEAN FIELD ANALYSIS 

In this appendix we review a mean field analysis for 
the q — 2 case following Refs. [s^. IsH . Ignoring vortex 
gauge fluctuations in Eq. (22), we can write down the 
vortex action with the effective magnetic field S — 1/2 in 
a lattice version 

Sm= 1 dr[^|9.$„p_^$tij;^ci>,„ + y(|$|)], 

n nm 

(A.l) 

where n, m label sites in the dual lattice, and the sign of 
the hopping integral t around a plaquette is — 1 owing 
to the TT fiux background. 

The vortex hopping term can be easily diagonalized 
in the eigenvectors Xn — (1 + V^) — e*^"" and Xn = 
e""^ [(l-|-\/2) + e"^"''], resulting in two low energy modes 
near the momentum k = (0,0) and k = (7r,0). Then, 
the low energy dynamics of this system can be described 
by the low energy vortex fields and ^Pin in <!>„ — 

In order to construct the effective vortex potential one 
can introduce the following two complex fields 

V-o^^o + i*!, Vi = *o-i*i- (A.2) 
Then, the symmetry transformations are given by 

Tx-.-ipo — >-ipi, i>i — >'ipo, 

Ty : ■00 — > #1, ■01 — > -#0, 

R./2 ■■ 00 e"/Vo, i'l — > e-"/Vi,(A.3) 

where T^^y) and are associated with lattice transla- 
tions and rotations. The LGW effective potential allowed 
by these symmetry operations is obtained to be 34] 

+«4|VonV'i|'-^^8[(V'o^i)' + i?-c.], (A.4) 

where is an effective vortex mass, U4 a local interac- 
tion, W4 the cubic anisotropy, and ug breaking the U(l) 
phase transformation V'o(i) ~^ ^^'''"''''^ i^o{i) ■ 

One cautious reader may ask how the coefficients in 
the LGW free energy functional can be determined. Al- 
though the symmetry constraints restrict the functional 
form of the effective potential, they cannot determine the 
remaining parameters in the free energy. Our question is 
whether these parameters are doping dependent or not. 
Remember that there is no doping dependence in the 
original vortex mass in Eq. (15). It depends on only 
the parameter U/t. However, the effective parameters for 



11 



the low energy vortex fields should be considered to be 
doping dependent. Consider a vortex vacuum resulting 
from large effective magnetic fields 6 in spite of < 
{U/t > {U/t)c), corresponding to chargon condensation. 
Decreasing hole concentration, the flavor number q of low 
energy vortices would increase. Decreasing hole concen- 
tration further, some components of the low energy vor- 
tices are expected to be condensed. In this respect the 
coefficients in the LGW free energy of (or ijji) can be 
considered to depend on hole concentration effectively. 

Based on the effective vortex potential Eq. (A4), one 
can perform a mean field analysis. Condensation of vor- 
tices occurs in the case of < and U4 > 0. The 



signs of V4 and vs then determine the ground state. For 
V4 < 0, both vortices have a nonzero vacuum expecta- 
tion value \ — \{ipi) \ 7^ 0, and their relative phase is 
determined by the sign of vg- In the case of > the 
resulting vortex state corresponds to the columnar dimer 
order, breaking the rotational and translational symme- 
tries. In the case of < the resulting phase exhibits 
the plaquette pattern, braking the rotational symmetries. 
On the other hand, if V4 > 0, the ground states are given 
by either |(^o)| ^ 0,J(i^i)| = or |(i^o)| = 0, KV'i)! ^ 0, 
and the sign of vs is irrelevant. In this case an ordi- 
nary charge density wave order at wave vector (tt, tt) is 
obtained, breaking the translational symmetries. 
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